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基于 无 人 机 高 光谱 遥感 的 烤烟 叶片 叶绿素 含量 估 测 


MEK SU, Æ A, 丰 , EZ, EO oh, KF, Tee 


(河南 农业 大 学 烟草 学 院 / 国 家 烟草 栽培 生理 生化 研究 基地 /烟草 行业 烟草 栽培 重点 实验 室 ,河南 郑州 450002) 


摘 要 : [目的 /意义 ] 烤烟 叶片 叶绿素 含量 (Leaf Chlorophyll Content, LCC) 是 表征 烤烟 光合 作用 、 营 养 状况 和 
长 势 的 重要 指标 。 本 研究 的 目的 为 高 效 精确 地 估 测 不 同 生 长 期 烤烟 LCC。[ 方 法 ] 以 中 烟 100 烟 叶 为 研究 对 象 ， 利 
用 无 人 机 搭载 Resonon Pika 工 高 光谱 成 像 仪 采 集 烤 烟 在 6 个 关键 生育 期 冠 层 反射 率 数据 。 基 于 相关 分 析 筛 选 了 21 种 
LCC 的 敏感 光谱 指数 ， 通 过 比较 不 同 光谱 组 合 及 不 同 回归 分 析 算 法 的 预测 精度 ， 最 终 建立 了 基于 多 种 光谱 指数 组 
AY LCC 回归 舍 测 模型 。 采 用 一 元 线性 回归 (Unary Linear Regression，ULR)、 多 元 线性 回归 (Multivariable Linear 
Regression, MLR), 、 偏 最 小 二 乘 回归 (Partial Least Squares Regression, PLSR) 、 支 持 向 量 回归 (Support Vector Re- 
gression, SVR) 和 随机 森林 回归 (Random Forest Regression, RFR) 5 种 建 模 方法 进行 LCC 估 测 。[ 结 果 和 讨论 ] 
在 不 同 生 育 期 大 部 分 光谱 参数 与 LCC 的 相关 性 达到 极 显著 (P<0.01); 相 较 于 传统 植被 指数 ， 新 组 合 的 光谱 指数 显 
著 提 升 了 与 LCC 的 相关 性 ; 对 单 变量 LCC 佑 测 模型 ULR， 以 移 栽 后 75 d 新 组 合 的 归 一 化 光谱 指数 与 红 光 比率 光谱 
指数 的 单 变 量 建 模 精度 最 高 ， 两 者 决定 系数 (Coefficient of Determination, R) 和 均 方 根 误差 (Root Mean Square 
Error, RMSE) 分 别 为 0.822 和 0.814，0.226 和 0.230。MLR、PLSR、SVR 和 RFR 建 模 方法 预测 结果 表明 RFR 
法 在 LCC 佑 测 中 效果 最 好 ， 甚 中 使 用 移 栽 后 75 d 数 据 验 证 集 的 尾 和 RMSE 可 达 0.919 和 0.146。[ 结 论 ] 本 研究 通过 
分 析 多 种 光谱 指数 与 烤烟 LCC 的 响应 规律 ， 构 建 可 靠 的 烤烟 叶片 LCC 估 测 模型 ， 可 为 烤烟 叶 LCC 佑 测 以 及 烤烟 的 
生长 发 育 监测 提供 理论 依据 和 技术 支撑 。 
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1 引 &8 费力 ， 效 率 低 下 。 相 较 于 采用 卫星 天 感 和 近 地 面 高 
光谱 人 工 采 集 方式 ， 无 人 机 遥感 具有 分 辩 率 高 、 效 


叶绿素 是 作物 叶片 中 主要 的 光合 色素 ， 叶 片 中 
叶绿素 含量 (Leaf Chlorophyll Content，LCC) 是 see NA UR ie ANGE AR (ee eT 
RUE TERRE A AE EULER AAAI PGES MA. 
标 。 烟 草 作为 重要 经 济 作物 ， 以 叶片 为 主要 收获 对 研究 表明 ， 基 于 植被 指数 《Vegetation Index, 
象 ， 因 此 叶绿素 含量 监测 非常 重要 。 传 统 的 烟草 农 VD 的 遥感 监测 具有 时 空 连续 性 ， 方 便 对 作物 长 势 
业 监 测 受 田 间 管 理 模式 及 生产 设备 的 限制 ， 费 时 、 ”信息 进行 长 期 监测 ， 在 叶绿素 估 测 中 具有 较 高 的 应 


收 稿 日 期 2023-03-14 

基金 项 目 : 烟草 行业 烟草 栽培 重点 实验 室 项 目 (30800665); 河南 省 科技 攻关 项 目 (172102110168) 
作者 简介 : 赖 佳 政 ,研究 方向 为 烟草 信息 学 。 E-mail: laijiazheng23@163.com 

“通信 作者 : 1 Æ, 博士， 讲师 ， 研 究 方向 为 烟草 信息 学 。E-mail: Zhangqian225@henau.edu.cn; 
2. 叶 协 锋 ， 博 士 ， 教 授 ， 研 究 方向 为 烟草 栽培 生理 和 健康 土壤 培育 。E-mail: yexiefeng@163.com 


ChinaXiv 合 作 期 刊 


69 


202308.00173v1 


ChinaXiv 


用 价值 "。 冯 海 宽 等 ”通过 分 析 VI 与 小 麦 叶 绿 素 
相对 含量 (Soil and Plant Analyzer Development, 
SPAD) 的 相关 性 ， 将 筛选 的 线性 组 合 指数 、 植 物 
生化 指数 、 归 一 化 指数 以 及 绿 光 归 一 化 指数 作为 偏 
最 小 二 乘 回归 (Partial Least Squares Regression, 
PLSR) 模型 的 输入 变量 ， 与 冬小麦 SPAD 值 联 合 构 
建 SPAD-PLSR 估 测 模型 ， 最 终 确 定 开 花期 的 SPAD 
估 测 精度 最 好 建 模 集 和 验证 集 民 分 别 为 0.78 和 
0.81。Qiao 等 ”研究 在 不 同 覆 盖 度 情况 下 大 田 玉米 
的 LCC 和 植被 指数 的 响应 特征 ， 采 用 随机 森林 
(Random Forest, RF) 和 偏 最 小 二 乘 两 种 方法 ， 构 
建 了 冠 层 LCC 佑 测 模型 。 田 明 璐 等 “选取 了 28 种 
与 作物 叶绿素 高 相关 的 植被 指数 ， 对 比 发 现 由 多 个 
光谱 参数 构建 的 PLSR 模型 对 棉花 LCC 佑 测 效果 最 
佳 。 研 究 人 员 尝 试 将 高 光谱 数据 对 应 的 波段 两 两 随 
机 组 合 ， 形 成 新 的 组 合 光谱 指数 用 以 反 演 作物 长 势 
fim. PME 采用 任意 两 波段 组 合 的 光谱 指数 
构建 了 马 铃 莫 LCC 估 测 模型 ， 相 较 于 其 他 传统 植被 
指数 ， 由 586 和 498 nm 组 合 的 归 一 化 光谱 指数 以 及 
586 和 462 nm 组 合 的 比率 光谱 指数 构建 的 反 演 模型 
精度 和 稳定 性 更 为 出 色 。Jia 等 “ 的 研究 表明 ， 新 
组 合 的 光谱 指数 有 助 于 敏感 波段 的 确定 ， 增 强 叶片 
须 含 量 的 相关 性 ， 在 烤烟 叶片 须 含 量 回归 模型 的 性 
能 表现 也 是 可 靠 的 。 近 年 来 许多 学 者 结合 机 器 学 习 
算法 预测 作物 农学 参数 ， 模 型 精度 和 重 棒 性 取得 较 
为 理想 的 表现 。 陈 晓 凯 等 "通过 不 同 光谱 变换 下 
的 最 优 光 谱 指 数 结合 BP 神经 网 络 和 RF 构建 冬小麦 
叶 面 积 指数 的 预测 模型 ，RF 预测 结果 优 于 BP 神经 
网 络 ， 并 且 很 好 地 避免 佑 测 模型 出 现 过 度 拟 合 和 拟 
合 不 足 的 情况 。 王 玉 娜 等 ”采用 随机 和 森林 回归 
(Random Forest Regression, RFR) 手段 估计 了 冬 
小 麦 的 损人 营养 指数 ， 得 到 了 性 能 理想 的 模型 (R= 
0.79, RMSE-0.13), Yang 等 "研究 结果 显示 采用 
RFR 可 以 显著 提高 马铃薯 块 蕉 形成 期 和 块 蕉 膨大 期 
的 地 上 部 生物 量 估计 精度 ， 同 时 显著 降低 光谱 数据 
的 共 线性 和 宛 余 度 。 

目前 已 发 表 无 人 机 遥感 相关 研究 多 是 针对 玉 
米 、 冬 小 麦 等 大 田 作 物 ， 有 关 烤 烟 相 关 人 研究 较 少 。 
通过 高 光谱 数据 进行 烟草 冠 层 LCC 的 定量 分 析 具 有 
重要 意义 和 应 用 价值 。 为 实现 烤烟 LCC 的 快速 估 


测 ， 本 研究 利用 无 人 机 获取 6 个 不 同 移 栽 生育 期 烤 
烟叶 高 光谱 影像 ， 通 过 筛选 光谱 参数 ， 利 用 一 元 线 
性 回归 (Unary Linear Regression, ULR), ZJC% 
性 回归 (Multivariable Linear Regression, MLR) 和 
PLSR = FRE St BPI, SCF) AAA (Support 
Vector Regression, SVR) 和 RFR 两 种 机 器 学 习 建 
模 方 法 分 别 构建 烤烟 的 LCC 估 测 模型 ， 并 获得 最 优 
佑 测 模型 ， 为 利用 无 人 机 遥感 开展 田间 烤烟 长 势 信 
息 快 速 、 无 损 监测 提供 技术 支持 。 


2 材料 与 方法 


2.1 试验 材料 


试验 田 位 于 河南 省 平顶山 市 部 县 植 烟 区 
(112*14'—113?45'E, 33?08'—34920'N), ， 属 于 温带 
大 陆 性 季风 和 气候， 气候 温 和 ， 雨 量 充 沛 ， 年 日 照 率 
53%， 平 均 气 温 15.7"C ， 无 霜 期 220 天 ， 常 年 降水 
量 678.5 mm. 

烤烟 品种 为 中 烟 100。 于 2022 年 5 月 8 日 移 栽 ， 
所 有 肥料 均 作 基肥 一 次 性 施 入 ,不 再 追肥 。 设 置 
5 个 水 平 氮 肥 梯 度 ，N0~N4 处理 分 别 为 0.0、24.0、 
48.0、72.0、96.0 kg/hm2， 每 个 处 理 设置 三 次 重复 ， 
共 15 个 小 区 ， 每 个 小 区 面积 约 为 144m2 (15.2 mX 
9.5 m), ， 如 图 1 所 示 。 其 他 栽培 和 管理 方法 按照 当 
地 优质 烟草 栽培 方式 进行 。 


2.2 数据 获取 


2.2.1 无 人 机 高 光谱 影像 获取 

获取 烤烟 6 个 生育 期 32、48、61、75、89 和 
109 d 的 影像 。 选 择 中 午 睛 朗 、 无 风 、 无 云 的 天 气 
条 件 ， 采 用 大 疆 M600 Pro 7 WERE IC AHL, RAR 
荷 5kg， 最 大 起 飞 重量 15.5 kg， 水 平 状态 下 最 大 的 
飞行 速度 能 达到 18 m/s， 续 航 时 间 约 30 min。 搭 载 
Resonon Pika 工 高 光谱 成 像 仪 ， 包 含 150 个 波段 的 
光谱 数据 ， 主 要 参数 如 表 1 所 示 。 地 面试 验 样品 采 
集 与 无 人 机 飞机 作业 同步 进行 。 

试验 于 上 述 烤 烟 生 育 期 11 : 00—14 : 00 进 行 测 
定 。 无 人 机 飞行 作业 时 ， 需 经 过 布置 在 试验 地 的 靶 
标 布 便于 后 期 辐射 定 标 处 理 数据 ,设置 航 高 100 m, 
飞行 速度 3 m/s， 最 小 定时 拍照 间隔 1.0 s， 航 向 重 
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注 : 该 图 由 自然 资源 部 标准 地 图 服务 网 站 下 载 的 审 图 号 为 GS(2019)1822 号 标准 地 图 制作 , 底 图 无 修改 
图 1 河南 省 平顶山 市 烤烟 试验 地 研究 区 分 布 图 


Fig. 1 Distribution map of the tobacco experimental study area in Pingdingshan city, Henan province 


#1 Pika 上 高 光谱 成 像 仪 主要 参数 
Table 1 Main parameters of Pika L hyperspectral imager 


649 nm 处 的 吸光 度 值 计 算 "“。LCC 计 算 方法 如 公 
式 (1) ~ (3) 所 示 。 


参数 名 称 数值 Ca = 13.95 X A4—76.88 X A, (1) 
光谱 范围 /nm 400~1000 Cb = 24.96 X A ,71.32 X A ges (2) 
光谱 分 辩 率 /nm 23 Ct = (Ca*Cb) X V/1000 X W (3) 
采样 间隔 /nm 1.07 其 中 ，Ca 和 Cb 分 别 表示 叶绿素 a 浓度 和 叶 
光谱 通道 数 /个 561 RA DIRE, mg/L; VRR ERRER, mL 
空间 通道 数 /个 900 下 表示 叶片 的 质量 ，g; As 和 As 分 别 为 665 和 
d hcic 249 649 nm 处 的 吸光 度 值 。 
尺寸 /em X em X em 10.0X 12.5X 5.3 


RA S096, 55H] HE 70% , 
17mm, 
2.2.2 ”烤烟 LCC 的 实验 室 测 定 
分 别 获取 移 栽 后 32、48、61、75、89 和 109d 
共 6 个 生育 期 烤烟 冠 层 叶片 进行 LCC 测 定 。 每 个 小 
区 进行 随机 破坏 性 采样 3 株 ， 每 次 数据 采集 共 获取 
45 个 地 面 样本 。 
将 烤烟 冠 层 叶片 转移 到 实验 室 测定 LCC。 用 
2 mm 打 孔 器 收集 新 鲜 样品 ， 混 匀 后 称 量 0.1 g ET, 
再 浸入 15 ml 的 95% 乙醇 溶液 中 。 在 浸 提 液 制 备 完 
成 后 ， 将 其 保存 在 低温 、 避 光 和 密闭 的 容器 中 ,在 
黑暗 中 提取 叶绿素 24 h。 经 过 上 暗 处 理 后 ， 叶 子叶 
白 绿 色 ， 使 用 Spectrumlab 722PC 比 色 分 光 光 度 计 
测量 叶片 色素 密度 。 提 取 的 LCC (Ct) 由 665 和 


镜头 垂直 朝 下 ， 焦 距 


2.3 光谱 参数 的 选 定 


为 了 降低 土壤 背景 、 大 气 水 分 等 环境 因素 对 冠 
a ， 大 多 数 研 究 人 员 通 过 利用 VI 增 
强 作物 光谱 信息 的 特征 并 且 实 现 高 光谱 数据 信息 降 
Hg em 本 研究 将 通过 寺 两 方面 选取 构建 烤烟 LCC 
佑 测 模型 的 光谱 指数 : 一 是 基于 已 发 表 的 对 作物 
LCC 敏感 的 18 种 植被 指数 ; 二 是 基于 波段 400~ 
1000 nm 范围 内 任意 两 个 波段 随机 组 合 的 差 值 光谱 
指数 (Difference Spectral Index, DSI), 、 比 值 光 谱 
指数 (Ratio Spectral Index, RSI) 和 归 一 化 光谱 指 
数 (Normalized Spectral Index, NDSI), ， 同 时 计算 
3 普 指数 与 LCC 的 相关 性 并 绘制 相关 性 等 势 

， 最 终 筛 选 相 关系 数 最 优 的 波段 组 合 为 组 合 光 谱 
pa 进行 研究 。 所 选取 的 光谱 参数 如 表 2 
所 示 。 
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表 2 光谱 参数 和 计算 公式 
Table 2 Spectral indices and formula of calculation 
光谱 参数 计算 公式 参考 文献 
植物 色素 比率 指数 (Plant Pigment Ratio ,PPR ) PPR=(Rss. -Rs0)/(Rssot Ryso) (4) [14] 
红 边 归 一 化 指数 (Red edge Normalized Difference Vegetation Index,RNDVI) = RNDVIF(R,,,—R.5)/( Rast Rigs) (5) [15] 
红 光 比率 光谱 指数 (Red Light Ratio Spectral Index, NIR) NIR-R,R, yy (6) [16] 
归 一 化 叶绿素 指数 (Normalized Difference Chlorophyll Index NDCI) NDCI-(R, -R,, )/( R4 R,,) (7) [17] 
简单 比值 指数 (Simple Ratio Index, SR) SR-R,,/R.,. (8) [15] 
红 光 叶绿素 光谱 指数 (Red-edge model index , CI, cage) Chien Reor Roa) 1 (9) [18] 
优化 土壤 调节 植被 指数 (Optimal Soil Adjusted Vegetation Index, OSAVI) OSAVIF1.16(Rggp—Regy )/( Rego + Royo t0- 16) 10) [19] 
红 边 归 一 化 光谱 指数 (Normalized Difference Red Edge Index, NDRE) NDRE-(R,,,-R,,)/( R4, FR, ) 11) [20] 
修正 归 一 化 光谱 指数 (Modified Normalized Difference Spectral Index, mND,,.) mND,,,7(R,, - Rj, )/(R;, FR, '42R,,) 12) [21] 
Vogelmann £1.11 48 4 ( Vogelman Red Edge Index, VOG) VOGER 46/R5 9 13) [22] 
红 边 位 置 指数 (Red Edge Position Index , REP) REP=700+40( (RootR;s0)/2-Ry00)/ (RRyo) 14) [23 
S 归 一 化 光谱 指数 550(ND,,,) ND, =R, — Resop)/ (Reso +R ssn) 15) [24] 
co 绿 光 叶绿素 光谱 指数 (Green Model Index, Cl...) CT een (Ryoo/Rsgo) 1 16) [18] 
p 修正 简单 比值 指数 (Modified Simple Ratio Index, mSRI) mSRI-(R,, -Rys )/(Rj,-R,,) 17) [25] 
o 转换 叶绿素 吸收 反射 率 指 数 (Transformed Chlorophyll Absorption in Reflec- TCARI-ACR, -R. J-02(R -R. XR.JR- G8) [26] 
e tance Index, TCARD 7 E dd C ID mE 
" 
TCARIOSAVI-[3((R,,-R,,)-0.20R, -R4 (Ra — 
ES Pena Rer) ) J/L 1.166 "T RS 16) ] E (19 ) si 
co 绿 光 归 一 化 植被 指数 (Green Normalized Difference Vegetation Index, GNDVI) GNDVI-(R,, -R,. )/( RR.) 20) [27] 
= 陆地 叶绿素 指数 (Meris Terrestrial Chlorophyll Index, MTCI) MTCI-(R,,,-R,,,)/( Raggi) 21) [28] 
N 差 值 植 被 指数 (Difference Vegetation Index, DSI) DSI-R.-R, 22) X 
> 比值 植被 指数 (Ratio Spectral Index , RSI) RSI=R,/R, 23) AS 
S< 归 一 化 植被 指数 (Normalized Difference Spectral Index, NDSI) NDSI=(R-R,)/(RAR,) 24) AM 
人 TE: i 和 j 表示 400~1000 nm 范围 内 的 任意 波段 位 置 ;， R AUR, 表示 在 上 和 7 位 置 的 反射 率 
= 2.4 模型 的 构建 与 检验 3 ”结果 与 分 析 
O 为 了 在 保证 模型 训练 充分 的 同时 ， 尽 可 能 减少 


模型 过 拟 合 的 风险 ， 以 建 模 集 和 验证 集 样本 比例 
2 : 1、 每 个 生育 期 45 个 样本 根据 比例 随机 划分 烤 
烟 样 本 数据 。 首 先 对 选 定 的 光谱 参数 和 LCC 进行 相 
关 性 分 析 以 筛选 对 LCC 敏感 的 光谱 指数 ， 然 后 将 得 
选 出 的 光谱 参数 作为 自 变 量 。 参 考 其 他 无 人 机 高 光 
谱 佑 测 作物 长 势 信息 所 选择 的 模型 及 常见 机 器 学 习 
模型 的 优势 ， 采 用 一 元 线性 回归 (Unary Linear Re- 
gression, ULR), 、MLR 和 了 PLSR 这 三 种 传统 线性 模 
型 ， 以 及 SVR 和 RFR 两 种 机 带 学 习 方法 对 单个 生 
育 期 和 全 生育 期 烤烟 LCC 进行 预 估 。 

采用 决定 系数 (Coefficient of Determination, 
RO) 和 均 方 根 误差 (Root Mean Square Error, 
RMSE) 对 模型 进行 精确 性 和 和 鲁 棱 性 检验 。 


3.1 烤烟 LCC 变化 


表 3 描 述 的 是 烤烟 冠 层 LCC 在 不 同 生育 期 的 变 
化 情况 (实验 室 测定 结果 )。 结 果 表 明 ， 烤 烟 生长 
前 期 ， 冠 层 LCC 会 迅速 增加 ， 在 团 棵 期 或 旺 长 期 
LCC 达 到 顶峰 ; 烤烟 进入 成 熟 期 后， 冠 层 叶 片 氮 素 
发 生 转 移 ，LCC 开始 减少 ， 而 烤烟 冠 层 LCC 在 现 
蕾 期 会 进一步 降低 。 从 整体 来 看 ， 冠 层 LCC 符合 先 
增 后 减 的 变化 趋势 ， 在 移 栽 后 75 d〈 旺 长 期 ) LCC 
达到 最 大 ， 之 后 开始 逐渐 下 降 。 在 各 生育 期 中 移 栽 
后 89 qd 变异 系数 最 大 ， 为 28.43%， 其 余 不 同 生育 期 
的 变异 系数 相差 不 大 ， 范 围 为 20.38% 一 25.26%， 各 
个 生育 期 的 变异 系数 均 低 于 全 生育 期 。 

全 生育 期 样本 范围 为 0.52~2.95 mg/g， 其 标准 
差 和 变异 系数 的 数值 表明 LCC 的 离散 程度 较 大 ， 反 
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映 出 试验 地 不 同 处 理 肥 力 之 间 是 有 差异 的 ， 也 保证 
了 佑 测 模型 在 一 定 范围 内 适用 。 

表 3 不 同 生育 期 烤烟 冠 层 叶片 LCC 的 变化 
Table 3 Variation of leaf chlorophyll content (LCC) in flue- 


cured tobacco leaves at different growth stages 


me Ts LCC Br HORR PEs Seem 
数 /个 (mg.g') f[i/(mgg?^) 数 /% 

移 栽 后 32d 45 0.73~2.19 1.53 0.35 23.20 
移 栽 后 48d — 45 0.84~2.23 1.61 0.33 20.38 
移 栽 后 61d 45 1.03~2.64 1.80 0.40 21.97 
移 栽 后 75d 45 0.87~2.95 2.05 0.52 25.26 
移 栽 后 89d | 45 0.79~2.46 1.60 0.45 28.43 
移 栽 后 109d | 45 0.52~1.49 0.88 0.20 22.13 
全 生育 期 270 | 0.52—2.95 1.58 0.52 33.22 


zr 


E. 全 生育 期 数据 集 包含 6 个 生育 期 所 有 样本 ;变异 系数 = 标准 
差 平 均值 


3.2 烤烟 LCC 和 光谱 参数 的 相关 性 分 析 


3.2.1 植被 指数 与 烤烟 LCC 的 相关 性 
将 6 个 生育 期 的 烤烟 LCC 与 植被 指数 进行 相关 
性 分 析 ， 结 果 如 表 4 所 示 。 


由 表 4 可 知 ， 除 移 栽 后 109 d， 其 余生 育 期 大 部 
分 植被 指数 和 LCC 呈 极 显著 相关 (p<0.01) 是 大 
部 分 植被 指数 与 叶绿素 是 正 相 关 。 其 中 ， 移 栽 后 
32 d4， 除 TCARI 和 了 PPR 表 现 出 不 相关 ， 其 余 均 为 极 
显著 相关 ， 最 大 相关 系数 是 0.771， 为 光谱 指数 
NDRE; 移 栽 后 48 d, TCARI 和 PPR 为 无 显著 相 
关 ，TCARIOSAVI 呈 显著 相关 (p<0.05) ， 其 余 为 
极 显 著 相 关 ，mND 作 为 相关 性 最 优 的 植被 指数 ， 
其 相关 系数 为 0.779; 移 栽 后 61 d 至 移 栽 后 89 d, 
这 三 个 生育 期 的 植被 指数 均 为 极 显 著 相 关 ，PPR、 
TCARI 和 TCARIOSAVI 表 现 为 极 显 著 负 相关 ， 但 
三 个 生育 期 相关 性 最 高 的 村 被 指数 各 不 相同 ， 移 栽 
后 61 d 相 关 性 最 高 为 mND,,， 移 栽 后 75 d 相 关 性 最 
高 为 NIR， 移 栽 后 89 d 相 关 性 最 高 为 TCARIOSA- 
VI， 其 相关 系数 分 别 为 0.870、0.902 和 0.877; BA 
后 109 d4， 除 了 PPR 、TCARI 和 TCARIOSAVI 表 现 
为 极 显著 负 相 关 ，REP 时 显著 相关 ， 其 余 植 被 指数 
无 显著 相关 。 全 生育 期 样本 中 所 有 植被 指数 和 LCC 
极 显著 相关 ，VOG 的 相关 性 最 大 ， 其 相关 系数 为 
0.771。 此 外 ， 在 各 个 生育 期 红 边 位 置 的 植被 指数 


表 4 不 同 生育 期 烤烟 LCC 与 植被 指数 的 相关 系数 


Table 4 Correlation coefficients between leaf chlorophyll content( LCC) in flue-cured tobacco and vegetation index at different 


growth stages 


植被 指数 移 栽 后 32 d 移 栽 后 48 d 移 栽 后 61 d 移 栽 后 75 d 移 栽 后 89 d 移 栽 后 109 d 全 生育 期 
NIR 0.587" 0.525" 0.698" 0.902" 0.822" 0.250 0.415" 
GNDVI 0.672" 0.7137 0.798" 0.848" 0.854" 0.184 0.426" 
NDRE 0.7717 0.675" 0.821" 0.883" 0.765" 0.249 0.7317 
ND, 5, 0.623" 0.734" 0.783" 0.8317 0.852" 0.179 0.426" 
mSRI 0.677" 0.744" 0.768" 0.835" 0.808" 0.229 0.708" 
PPR 0.068 —0.032 -0.565" -0.616" -0.7577 -0.624" -0.323" 
REP 0.587" 0.530" 0.746" 0.876" 0.849" 0.306 0.664" 
SR 0.698" 0.782" 0.858" 0.852" 0.808" 0.183 0.601" 
RNDVI 0.698" 0.792" 0.864" 0.8407 0.8277 0.212 0.530" 
VOG 0.738" 0.673" 0.775" 0.869" 0.694" 0.221 0.7717 
NDCI 0.482" 0.649" 0.697" 0.742" 0.847" 0.069 0.273" 
CI, 0.767" 0.674" 0.817” 0.896" 0.754" 0.219 0.729" 
CE 0.655" 0.725" 0.821" 0.864" 0.830" 0.157 0.490" 
mND,,, 0.685" 0.799" 0.870" 0.846" 0.821" 0.170 0.460" 
TCARI -0.249 0.087 -0.675" -0.709" -0.873" -0.4737 -0.620" 
TCARI/OSAVI -0.590" -0.311 -0.759" -0.7617 -0.8777 -0.404" -0.630" 
OSAVI 0.526" 0.751" 0.467” 0.387” 0.413" -0.262 0.161" 
MTCI 0.691" 0.657" 0.806" 0.866 0.8107 0.259 0.682" 


iE: * 表示 0.05 水 平 显 著 ，** 表 示 0.01 KEE. Fi 
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和 LCC 含 量 的 相关 性 均 优 于 其 他 植被 指数 ， 其 中 
NDRE, Clre, RNDVI, VOG 以 及 mND;,; 在 不 同 生 
育 期 的 相关 性 ， 除 了 移 栽 后 109 d， 相 关系 数 均 在 
0.67 以 上 ,， 均 呈 极 显著 相关 。 以 上 结果 表明 ，VI 对 
LCC 的 敏感 性 与 作物 生长 阶段 密切 相关 。 随 着 生育 
期 推进 ， 叶 片 中 的 叶绿素 会 发 生 转移 或 降解 ， 这 意 
味 着 有 些 VIs 不 能 用 于 估 测 LCC。 因 此 ， 确 定 不 同 
生长 阶段 VIs 对 LCC 的 敏感 性 对 于 基于 无 人 机 高 光 
谱 图 像 的 作物 参数 估计 非常 重要 。 

对 全 生育 期 而 言 ， 相 关系 数 绝对 值 较 大 的 植 
被 指数 里 绝 大 部 分 是 红 边 位 置 反射 率 的 植被 指数 。 
在 不 同 生育 期 相关 性 最 优 的 指数 分 别 为 NDRE 
(32d), mND,,, (48d, 61d), NIR (75d), TCA- 
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(a) f& RUG 32d 


RI/OSAVI (89d), TCARI (10d) 以 及 VOG (4 
生育 期 )。 分 析 原 因 可 能 是 因为 NDRE、mND;;、 
NIR 和 VOG 主要 使 用 红 边 位 置 和 近 红 外 波段 信息 ， 
这 类 植被 指数 对 于 烤烟 LCC 的 建 模 估 测 具有 较 好 的 
应 用 潜力 ， 参 与 建 模 可 以 有 效 提 升 模型 的 精确 度 和 
和 鲁 棒 性 ”1。 
3.2.2 组合 光谱 指数 与 烤烟 LCC 的 相关 性 

以 无 人 机 高 光谱 波段 范围 400~1000 nm 内 任意 
两 波段 组 合计 算得 到 的 光谱 指数 RSI、DSI 和 
NDSI， 并 与 叶片 叶绿素 进行 相关 性 分 析 。 计 算 并 
绘制 各 个 生育 期 和 全 生育 期 的 相关 性 等 势 图 ， 如 
图 2 和 图 3 所 示 。 
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(了 f) 移 栽 后 109 d 


图 2 不 同 生育 期 烤烟 LCC 和 组 合 光谱 指数 的 相关 性 等 势 图 


Fig. 2 Isopotential plots based on correlation between LCC in flue-cured tobacco at different growth stages and new combination of 


spectral indices 
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图 3 基于 全 生育 期 烤烟 LCC 和 组 合 光谱 指数 的 相关 性 等 势 图 
Fig. 3 Isopotential plots based on correlation between LCC in flue-cured tobacco leaves during the whole growth stages and new com- 


bination of spectral indices 
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结果 表明 ， 除 了 移 栽 后 109 d， 各 生育 期 的 敏 
感 波段 相对 较为 集中 ， 与 LCC 相 关 性 较 高 的 光谱 指 
数组 合 主 要 分 布 在 780~940 nm Fil 520—710 nm 之 
间 ， 移 栽 后 109 d 相 较 其 他 生育 期 效果 不 佳 。 三 种 
随机 组 合 的 光谱 指数 大 部 分 小 于 相关 系数 0.3801， 
达到 极 显 著 相 关 水 平 的 区 域 较 少 。 全 生育 期 的 敏感 
波段 相对 较为 分 散 ， 不同 光谱 指数 之 间 规 律 相 差 不 
大 ,但 相 较 于 传统 植被 指数 ， 其 相关 性 有 很 大 提 
升 。 以 光谱 指数 与 烤烟 LCC 相关 性 最 高 为 原则 ， 筛 
选 出 不 同 生 育 期 和 全 生育 期 相关 系数 最 优 的 21 种 
组 合 光 谱 指 数 ， 如 表 S$ 所 示 。 结 果 表 明 ， 相 较 于 已 
有 的 植被 指数 ， 组 合 光 谱 指 数 和 LCC 的 相关 性 有 着 


更 进一步 的 提升 ， 并 且 烤 烟 LCC 与 筛选 出 的 最 优 光 
谱 指 数 均 通过 0.01 水 平 的 显著 性 检验 。 对 比 三 种 组 
合 光 谱 指 数 ， 归 一 化 光谱 指数 与 烤烟 LCC 在 移 栽 后 
32、48、61、75 和 109 d 相 关 性 最 高 ， 这 些 生 育 期 
的 最 优 波段 组 合 为 797 和 719 nm、727 和 697 nm, 
736 和 706 nm, 775 和 745 nm 以 及 455 和 541 nm, 
其 相关 系数 分 别 为 0.795、0.805、0.878、0.912 和 
0.669; 组 合 差 值 光谱 指数 在 移 栽 后 89 d 时 与 LCC 
相关 性 最 高 ，DSI cases, ro) 与 烤烟 LCC 的 相关 系数 
为 0.892; 比值 光谱 指数 在 全 生育 期 数据 集中 相关 
系数 最 大 ，RSI qus sus 与 烤烟 LCC 的 相关 系数 为 
0.807. 


X5 烤烟 不 同 移 栽 生育 期 组 合 光 谱 指数 入 选 波段 和 LCC 相关 系数 


Table 5 Selected bands of new combination of spectral indices and LCC correlation coefficients at different growth stages 


移 栽 天 数 /d 光谱 指数 相关 系数 移 栽 天 数 /d 光谱 指数 相关 系数 
Si 0.754” DSI soo 0.892" 
32 RSL css p79) 0.792" 89 RSI pose. noy) 0.870" 
NDS ikana 0.795" ND ST aay 0.881" 
DSL iun 0.725" DSL gyioy 0.521" 
48 RSI ici) 0.7977 109 RSI. enu) 0.668" 
NDST 797.2697) 0.805" NDSI pass sui) 0.669" 
DSL s s 0.822" DSL iis 0.726" 
61 RIST a 0.878” 全 生育 期 RSI ci 0.807" 
NDST D kari) 0.878” ND Siati 0.806" 
DSL 5255 0.835" 
75 RO ea 0.912™ 
NDSI 0.912™ 


(R775,R745) 


3.3 烤烟 LCC 估 测 模型 的 构建 


3.3.1 ”基于 单个 自 变 量 烤烟 LCC 估 测 模型 

在 叶片 叶绿素 和 光谱 参数 相关 性 分 析 的 基 
础 上 上， 筛选 出 6 个 生育 期 和 全 生育 期 相关 系数 
最 高 的 植被 指数 NDRE、MND,,,、NIR、TCARI 
OSAVI, PPR 和 VOG， 光 谱 指 数 NDSI (o ro) > 
NDSI ro ~ NDSI as. ro ~ NDSI ies. es 


后 75 d， 其 建 模 集 和 验证 集 的 精度 是 逐渐 增加 的 ， 
移 栽 后 75 d 的 NIR 作 为 各 个 生育 期 佑 测 效 果 最 佳 的 
植被 指数 模型 ， 其 建 模 集 和 了 验证 集 的 民 分 别 是 
0.814 和 0.829，RMSE 分 别 是 0.230 和 0.253， 而 移 
栽 后 75 d 的 NDSI ans, sus 是 各 个 生育 期 估 测 精度 最 
佳 的 组 合 光 谱 指 数 ， 其 建 模 集 和 验证 集 的 RR 分别 是 
0.822 和 0.862, RMSE 分 别 是 0.226 和 0.227; 4E 
后 75 d BAB Bet 109 d 建 模 集 和 验证 集 的 精度 开 


DSI iess. rs) > RSI cass, xus 和 RSI ws， 分 别 与 
烤烟 叶片 叶绿素 建立 ULR 模型 ， 建 模 集 和 验证 集 
的 估 测 结果 如 表 6 所 示 。 各 个 生育 期 的 组 合 光 谱 指 
数 无 论 是 建 模 集 还 是 验证 集 ， 其 估 测 精度 大 多 优 于 
传统 植被 ， 但 不 同 生育 期 的 光谱 参数 的 表现 存在 差 
异 。 植 被 指数 和 优化 光谱 指数 从 移 栽 后 32 d 至 移 栽 


始 下 降 ， 而 在 移 栽 后 109 d 时 ， 植 被 指数 和 优化 
组 合 光 谱 指 数 在 各 个 生育 期 中 估 测 LCC 的 效果 都 
是 最 差 。 相 较 于 其 他 单个 生育 期 ， 由 不 同 生育 期 
270 个 样本 组 成 的 全 生育 期 数据 集 估 测 表现 较为 不 
理想 ，VOG 建 模 集 和 验证 集 的 尼 分 别 是 0.602 和 
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0.558, RMSE 分 别 是 0.348 和 0.349, RSI arse. rm) 分 别 是 0.333 和 0.304。 
建 模 集 和 验证 集 的 尼 分 别 是 0.636 和 0.686，RMSE 
表 6 不 同 生 育 期 单个 光谱 参数 的 烤烟 LCC 估 测 模型 及 验证 


Table 6 Prediction models and validation of flue-cured tobacco LCC with single spectral parameter at different growth stages 
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建 模 集 验证 集 
移 栽 天 数 /d 光谱 参数 回归 方程 - - 
R RMSE Rt RMSE 
NDRE y710.528x-0.968 (25) 0.606 0.234 0.527 0.198 
= NDSI s y=10.912x-1.035 — (26) 0.617 0.231 0.662 0.167 
MG EU c c gn E E 
T NDSI 797.2697) y=9.709x-3.949 — (28) 0.646 0.189 0.738 0.163 
Ec TV. r 2 — EE MEM a —— ON 
ib NDSI. aie y=12.001x-3.651 — (30) 0.774 0.216 0.749 0.225 
s Qu C o 4 Er ER MEE ME 
i NDS kois Ris) y=63.199x-0.391 (32) 0.822 0.226 0.862 0.227 
ee 
DSI sp nady y=126.691x+4.004 (34) 0.772 0.214 0.847 0.169 
M eg RO MT EE M. E 
i RSI eer) y=7.509x-0.890 — (36) 0.473 0.131 0.462 0.149 
uu uso — 1 meme E EE E 
全 生育 期 
RSie) y=7.026x-7.873 — (38) 0.636 0.333 0.686 0.304 
3.3.2 ”基于 多 自 变量 烤烟 LCC 估 测 模型 RMSE 分 别 为 0.893、0.850 和 0.214、0.217。 


基于 高 光谱 遥感 的 作物 生理 生化 参数 估 测 的 研 
究 显 示 ， 相 较 于 构建 单一 光谱 参数 估 测 模型 ， 以 多 
个 光谱 参数 作为 自 变 量 的 建 模 方法 对 估 测 模型 精度 
和 和 鲁 棒 性 都 有 较为 明显 的 提升 707. 但 由 于 自 变 量 
之 间 存 在 共 线 性 的 问题 ， 在 构建 烧 烟 LCC 佑 测 模 型 
过 程 中 需要 选取 适当 数量 的 自 变 量 来 防止 预测 结 
过 度 拟 合 。 因 此 ， 本 研究 选择 每 个 生长 生育 期 以 及 
全 生育 期 数据 集中 相关 系数 最 高 的 7 个 植被 指数 和 
3 个 组 合 光 谱 指 数 作 为 MLR、PLSR、SVR 和 RFR 
四 种 建 模 方 法 的 自 变 量 。 各 个 生长 阶段 通过 不 同 的 
建 模 分 析 方 式 构建 烤烟 LCC 估 测 模型 ， 模 型 的 建 模 
集 和 验证 集 估 测 效果 如 表 7 所 示 。 

从 表 7 的 结果 来 看 ， 相 较 于 单个 光谱 参数 ， 多 
个 光谱 参数 的 模型 佑 测 精度 在 各 个 生长 阶段 和 全 
生育 数据 集中 有 着 大 幅度 的 提升 。 对 比 各 个 生育 
期 的 建 模 集 佑 测 精度 发 现 ， 在 移 栽 后 32、48、61 
和 109 d， 模 型 LCC-RFR 的 估 测 性 能 最 好 ，R 分 别 
为 0.733、0.744、0.866 和 0.678，RMSE 分 别 为 
0.186、0.165、0.176 和 0.120， 而 基于 LCC-MLR 的 
估 测 模型 ， 在 移 栽 后 75 和 89 d 效 果 最 优 ， 忆 和 


对 比 验证 集 估 测 精 度 来 看 ， 移 栽 后 32 d, PEA 
LCC-MLR 的 估 测 精度 最 佳 ，R A RMSE 分别 为 
0.679 和 0.215， 其 余 五 个 生育 期 模型 LCC-RFR AY fH 
WU ARE RE, RAIA 0.757, 0.832, 0.919, 0.842 
和 0.580, RMSE 分 别 为 0.173, 0.185, 0.146, 
0.246 和 0.131。 在 全 生育 数据 集中 LCC-RFR 模型 的 
估 测 精度 较 优 于 LCC-MLR、LCC-PLSR 和 LCC- 
SVR 模型 ， 对比 建 模 集 ，R 增加 了 19.06%、 
18.62% 和 29.51%，RMSE 降低 了 31.93%、29.51% 
和 28.24%; 对 比 验证 集 ，R? 增 加 了 8.21%、12.62% 
和 8.17%，RMSE 降 低 了 3.76%、9.33% 和 4.55%。 

综 上 所 述 ， 通 过 各 个 生育 期 的 建 模 集 和 验证 集 
的 佑 测 精度 对 比分 析 ，LCC-RFR 模型 的 估 测 效果 
最 为 理想 ，LCC-MLR 和 LCC-SVR 模 型 估 测 精度 次 
之 ，LCC-PLSR 模 型 估 测 精度 最 差 ; 不 同 生育 期 模 
型 预测 精度 和 稳定 性 由 高 到 低 分 别 为 移 栽 后 75、 
61、89、48、32 和 109 d。 

3.3.3 ”烤烟 LCC 空间 分 布 图 

根据 不 同 生育 期 构建 的 烤烟 LCC 佑 测 模 型 结 

可 知 ，LCC-RFR 佑 测 模 型 的 精度 和 稳定 性 相对 最 
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m7 不 同 生育 期 多 个 光谱 参数 的 烤烟 LCC 估 测 模型 
Table 7 LCC prediction models of flue-cured tobacco with mul- 


tiple spectral parameters at different growth stages 


建 模 集 验证 集 


R RMSE R RMSE 


LCC-MLR 0.673 0.254 0.679 0.215 
LCC-PLSR 0.590 0.247 0.626 0.162 


2 LCC-SVR 0.647 0.232 0.569 0.176 
LCC-RFR 0.733 0.186 0.660 0.210 
LCC-MLR 0.680 0.240 0.742 0.148 
LCC-PLSR 0.607 0.217 0.621 0.192 

si LCC-SVR 0.702 0.680 0.691 0.164 
LCC-RFR 0.744 0.165 0.757 0.173 
LCC-PLSR 0.746 0.238 0.750 0.227 

? LCC-SVR 0.830 0.197 0.736 0.228 
LCC-RFR 0.866 0.176 0.832 0.185 
LCC-PLSR 0.766 0.293 0.807 0.204 

? LCC-SVR 0.856 0.204 0.824 0.269 
LCC-RFR 0.891 0.205 0.919 0.146 
LCC-PLSR 0.765 0.221 0.799 0.224 

" LCC-SVR 0.790 0.218 0.810 0.184 
LCC-RFR 0.846 0.163 0.842 0.246 

T LCC-PLSR 0.484 0.119 0.429 0.192 


LCC-SVR 0.513 0.117 0.606 0.168 
LCC-RFR 0.678 0.120 0.580 0.131 
LCC-MLR 0.718 0.303 0.741 0.274 
LCC-PLSR 0.720 0.292 0.712 0.291 


LCC-SVR 0.731 0.287 0.742 0.276 
LCC-RFR 0.854 0.206 0.802 0.264 


佳 。 因 此 ， 本 人 研究 在 消除 土壤 背景 和 杂 草 的 干扰 
后 ， 利 用 LCC-RFR 估 测 模型 选择 烤烟 三 个 生育 期 
进行 无 人 机 高 光谱 影像 的 反 演 填 图 ， 得 到 了 移 栽 后 
48、61 和 75 d 的 烤烟 LCC 空间 分 布 图 ， 如 图 4 所 
示 。 整 体 来 看 ， 从 移 栽 后 48 到 75 d， 随 着 生育 期 的 
推进 颜色 逐渐 加 深 ,表明 烤烟 LCC 越 来 越 大 。 不 同 
重复 之 间 也 有 所 差异 ， 烤 烟 长 势 分 布 由 强 到 弱 为 
“B=” “He” BE". MPR BAS d, 
LCC WYSE EIE 14~1.9 mg/g 之 间 ， 从 移 栽 后 61 dJF 


始 ， 重 复 一 与 重复 二 、 重 复 三 LCC 的 差异 愈 发 明 
fh, 重复 一 LCC 估 测 的 范围 大 多 在 1.8~2.3 mg/g, 
重复 二 和 重复 三 估 测 LCC 大 多 在 2.1~2.5 mg/g it 
围 内 。 移 栽 后 75 d LCC 整体 分 布 较 密 ，LCC 的 范 
EE 2.2—3.1 mg/g 之 间 ， 这 可 能 与 烤烟 生长 特性 有 
关 ， 在 烤烟 生育 前 期 至 旺 长 期 ， 烤 烟 LCC 呈 逐渐 上 
升 趋势 。 通 过 估 测 多 个 生育 期 的 LCC， 在 无 人 机 高 
光谱 图 像 上 反 演 的 LCC 空间 分 布 与 实测 结果 相符 ， 
这 表明 无 人 机 搭载 高 光谱 成 像 仪 可 提供 小 范围 内 烤 
烟 LCC 的 空间 分 布 信息 ， 具 有 较 高 的 分 辨 率 和 精 
度 ， 可 以 较 好 地 分 辨 出 烤烟 各 个 生育 期 烤烟 长 势 


4 讨 论 


本 研究 通过 将 烤烟 LCC 与 选取 的 光谱 参数 进行 
相关 分 析 ， 结 果 表 明 本 研究 构建 的 组 合 光 谱 指 数 与 
LCC 的 相关 性 明显 高 于 传统 植被 指数 ， 且 基于 单个 
光谱 构建 LCC 回归 模型 中 ， 前 者 的 估 测 精度 和 稳定 
性 在 各 个 生育 期 包括 全 生育 期 都 优 于 后 者 。 前 人 研 
究 表明 构建 组 合 光谱 指数 能 够 有 效 减 少 外 界 环 境 因 
素 干扰 ， 精 准 识别 作物 生长 参数 的 敏感 波段 ， 对 农 
作物 气 素 、 叶 绿 素 和 叶 面 积 指数 的 监测 精度 以 及 模 
型 鲁 棱 性 有 显著 提升 “*"。 但 也 有 学 者 认为 ， 组 
合 光 谱 指 数 所 包含 的 光谱 信息 不 足 ， 不 能 充分 利用 
从 可 见 光 到 近 红 外 波段 的 高 光谱 数据 信息 T. A 
此 ， 本 研究 采用 传统 植被 指数 和 优化 光谱 指数 相 结 
合 的 方式 试图 构建 更 为 综合 的 烤烟 LCC 佑 测 模型 。 
试验 结果 表明 ， 基 于 多 个 光谱 参数 组 合 构建 的 烤烟 
LCC 佑 测 模型 精度 和 拟 合 性 都 优 于 单个 光谱 参数 构 
建 的 模型 。 多 个 自 变量 可 以 有 效 避 免 光 谱 信 息 包含 
不 足 的 情况 ， 充 分 利用 高 光谱 数据 ， 但 是 光谱 参数 
过 多 也 容易 出 现 共 线 性 和 过 拟 合 的 问题 ， 因 此 , 今 
后 研究 需要 关注 对 烤烟 LCC 估 测 模型 自 变 量 的 
控制 。 

此 外 ,不 同 生育 期 的 模型 精度 各 不 相同 ， 且 随 
着 生育 期 的 推进 ， 模 型 精度 呈现 先 增 后 减 的 变化 趋 
势 ， 这 可 能 是 由 于 烤烟 在 生长 发 育 的 不 同 阶段 ， 伴 
随 着 冠 层 结构 的 变化 、 生 物 量 的 积累 ， 以 及 在 生殖 
生长 阶段 冠 层 叶 片 氮 含量 的 分 配 ， 导 致 光谱 参数 的 
预测 性 能 在 不 同 阶段 发 生变 化 。 对 比分 析 MLR, 
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图 4 烤烟 不 同 生育 期 的 LCC 空间 分 布 图 
Fig. 4 LCC spatial distribution map of flue-cured tobacco 


plants at different growth stages 


PLSR, SVR 和 RFR 四 种 多 自 变 量 LCC 估 测 模型 ， 
基于 RFR 模型 估 测 烤烟 LCC 的 精度 和 稳定 性 效果 


ChinaXiv 合 作 期 刊 


最 好 ， 主 要 原因 是 RFR 建 模 时 不 易 出 现 过 拟 合 现 
象 ， 能 够 有 效 降 低 异常 值 和 噪声 的 干扰 “。 利 用 
多 生育 期 数据 集 构建 全 生育 期 模型 会 降低 模型 预测 
精度 ”， 一 方面 可 能 是 因为 冠 层 结构 和 土壤 背景 
情况 随 着 生育 期 的 推进 在 不 断 改变 ， 由 冠 层 光谱 对 
LCC 的 响应 规律 发 生变 化 造成 ， 另 一 方面 ， 可 能 是 
由 不 同 生育 期 冠 层 光谱 的 测定 误差 全 加 造成 。 此 
外 ， 目 前 本 人 研究 的 数据 集 涵盖 了 不 同 生育 阶段 和 毛 
素 水 平 下 的 样本 ， 但 构建 的 组 合 光谱 指数 尚未 在 不 
同 环境 和 年 份 中 得 到 验证 。 为 了 提高 模型 的 稳健 性 
和 泛 化 能 力 ， 需 要 进一步 测试 不 同 烤烟 品种 的 适 
用 性 。 


5 结 it 


本 研究 探索 驳 中 烟 区 烤烟 LCC 与 无 人 机 高 光谱 
影像 的 响应 规律 ， 通 过 多 种 光谱 指数 与 烤烟 LCC 的 
相关 分 析 ， 构 建 烤烟 叶片 叶绿素 佑 测 含量 的 ULR、 
MLR、PLSR、SVR 和 RFR 模 型 ， 得 到 以 下 结论 : 

(1) 所 涉 大 部 分 植被 指数 和 LCC 是 极 显著 相 
关 。 其 中 ， 移 栽 后 32 一 109 d 相 关系 数 最 大 的 植被 
指数 分 别 是 NDRE MND,,. mND,,. NIR, TCA- 
RIOSAVI 和 REP， 全 生育 期 数据 集中 所 有 植被 指 
数 与 烤烟 LCC 均 呈 极 显著 相关 ，VOG 的 相关 性 最 
大 (r=0.771)。 

(2) 根据 随机 两 两 波段 组 合 的 光谱 指数 与 烤烟 
LCC 含 量 相关 性 最 大 原则 ， 筛 选 出 不 同 生育 期 和 全 
生育 期 相关 系数 最 高 的 21 种 光谱 指数 ， 传 统 植被 
虽 数 和 组 合 光谱 指数 与 LCC 的 敏感 波段 大 多 集中 在 
红 边 位 置 ， 但 整体 上 组 合 光谱 指数 的 相关 性 较 前 者 
更 为 出 色 ， 其 中 移 坊 后 75 d 的 NDSI wuws 和 
RSI ans, vus 的 相关 性 最 高 ， 相 关系 数 r 均 为 0.912。 

(3) 在 叶片 叶绿素 和 光谱 参数 相关 性 分 析 的 基 
础 上 ， 筛 选 出 6 个 生育 期 和 全 生育 期 烤烟 LCC 相关 
性 最 高 的 光谱 参数 与 烤烟 叶片 叶绿素 分 别 建立 ULR 
模型 ， 不 同 移 栽 生育 期 的 优化 光谱 指数 无 论 是 建 模 
集 还 是 验证 集 其 估 测 精度 均 较 优 于 传统 植被 指数 。 

(4) 相 较 于 单个 光谱 参数 ， 基 于 多 个 光谱 参数 
构建 模型 的 估计 精度 在 各 个 生长 阶段 和 全 生育 数据 
集中 均 有 大 幅度 提升 ， 其 中 移 栽 后 75 d 由 LCC- 
RFR 构建 的 模型 最 佳 UER RA RMSE 可 达 


ChinaXiv (ERAT 


79 


202308.00173v1 


ChinaXiv 


0.919 #10.146, ARKE, RFR 模型 的 LCC 估 测 效 
果 最 为 理想 ，LCC-MLR 和 LCC-SVR 模型 估 测 精度 
次 之 ，LCC-PLSR 模型 估 测 精度 最 差 ; 不 同 生育 期 
模型 预测 精度 和 稳定 性 由 高 到 低 分 别 为 移 栽 后 75、 


61, 


89, 48, 3251109 d, 
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Monitoring of Leaf Chlorophyll Content in Flue-Cured 
Tobacco Based on Hyperspectral Remote Sensing of 
Unmanned Aerial Vehicle 
LAI Jiazheng, LI Beibei, CHENG Xiang, SUN Feng, CHENG Juting, 
WANG Jing, ZHANG Qian’, YE Xiefeng 
(College of Tobacco Science, Henan Agricultural University/ National Tobacco Cultivation and Physiology and Biochem- 
istry Research Center/ Key Laboratory for Tobacco Cultivation of Tobacco Industry, Zhengzhou 450002, China) 
Abstract: 


[Objective] Leaf chlorophyll content (LCC) of flue-cured Tobacco is an important indicator for characterizing the photosynthesis, nu- 


tritional status, and growth of the crop. Tobacco is an important economic crop with leaves as the main harvest object, it is crucial to 
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monitor its LCC. Hyperspectral data can be used for the rapid estimation of LCC in flue-cured tobacco leaves, making it of great sig- 
nificance and application value. The purpose of this study was to efficiently and accurately estimate the LCC of flue-cured tobacco 
during different growth stages. 

[Methods] Zhongyan 100 was chose as the research object, five nitrogen fertilization levels were set. In each plot, three plants were 
randomly and destructively sampled, resulting in a total of 45 ground samples for each data collection. After transplanting, the reflec- 
tance data of the flue-cured tobacco canopy at six growth stages (32, 48, 61, 75, 89, and 109 d ) were collected using a UAV equipped 
with a Resonon Pika L hyperspectral. Spectral indices for the LCC estimation model of flue-cured tobacco were screened in two ways: 
(1) based on 18 published vegetation indices sensitive to LCC of crop leaves; (2) based on random combinations of any two bands in 
the wavelength range of 400—1000 nm. The Difference Spectral Index (DSI), Ratio Spectral Index (RSI), and Normalized Spectral In- 
dex (NDSI) were calculated and plotted against LCC. The correlations between the three spectral indices and leaf LCC were calculat- 
ed and plotted using contour maps. Five regression models, unary linear regression (ULR), multivariable linear regression (MLR), par- 
tial least squares regression (PLSR), support vector regression (SVR), and random forest regression (RFR), were used to estimate the 
chlorophyll content. A regression estimate model of LCC based on various combinations of spectral indices was eventually construct- 
ed by comparing the prediction accuracies of single spectral index models multiple spectral index models at different growth stages. 
[Results and Discussions] The results showed that the LCC range for six growth stages was 0.52—2.95 mg/g. The standard deviation 
and coefficient of variation values demonstrated a high degree of dispersion in LCC, indicating differences in fertility between differ- 
ent treatments at the test site and ensuring the applicability of the estimation model within a certain range. Except for 109 d after 
transplanting, most vegetation indices were significantly correlated with LCC (p<0.01). Compared with traditional vegetation indices, 
the newly combined spectral indices significantly improved the correlation with LCC. The sensitive bands at each growth stage were 
relatively concentrated, and the spectral index combinations got high correlation with LCC were mainly distributed between 780— 
940 nm and 520—710 nm. The sensitive bands for the whole growth stages were relatively dispersed, and there was little difference in 
the position of sensitive band between different spectral indices. For the univariate LCC estimation model, the highest modeling accu- 
racy was achieved using the newly combined Normalized Spectral Index and Red Light Ratio Spectral Index at 75 d after transplant- 
ing. The coefficients of determination (R°) and root mean square errors (RMSE) for the modeling and validation sets were 0.822, 
0.814, and 0.226, 0.230, respectively. The prediction results of the five resgression models showed that the RFR algorithm based on 
multivariate data performed best in LCC estimation. The R’ and RMSE of the modeling set using data at 75 d after transplanting were 
0.891 and 0.205, while those of the validation set reached 0.919 and 0.146. In addition, the estimation performance of the univariate 
model based on the whole growth stages dataset was not ideal, with R’ of 0.636 and 0.686, and RMSE of 0.333 and 0.304 for the mod- 
eling and validation sets, respectively. However, the estimation accuracy of the model based on multiple spectral parameters was sig- 
nificantly improved in the whole growth stages dataset, with R^ of 0.854 and 0.802, and RMSE of 0.206 and 0.264 for the modeling 
and validation sets of the LCC-RFR model, respectively. In addition, in the whole growth stages dataset, the estimation accuracy of 
the LCC-RFR model was better than that of the LCC-MLR, LCC-PLSR, and LCC-SVR models. Compared with the modeling set, R° 
increased by 19.06%, 18.62%, and 29.51%, while RMSE decreased by 31.93%, 29.51%, and 28.24%. Compared with the validation 
set, R’ increased by 8.21%, 12.62%, and 8.17%, while RMSE decreased by 3.76%, 9.33%, and 4.55%. 

[Conclusions] The sensitivity of vegetation indices (VIs) to LCC is closely connected to the tobacco growth stage, according to the 
results this study, which examined the reaction patterns of several spectral indices to LCC in flue-cured tobacco. The sensitivity of VIs 
to LCC at various growth stages is critical for crop parameter assessment using UAV hyperspectral photography. Five estimation mod- 
els for LCC in flue-cured tobacco leaves were developed, with the LCC-RFR model demonstrating the greatest accuracy and stability. 
The RFR model is less prone to overfitting and can efficiently decrease outlier and noise interference. This work could provide theoret- 


ical and technological references for LCC estimate and flue-cured tobacco growth monitoring. 
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